Dephasing in quantum dot molecules via exciton-acoustic phonon coupling 
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We develop a theory of the linear optical spectrum of excitons in quantum dot molecules, including 
the effect of exciton-phonon coupling beyond the Markov limit. The model reproduces the general 
trend of the zero-phonon line broadening as a function of interdot distance, that were recently 
measured. The unexpectedly broad linewidths and their large variation are explained in terms of 
both the non-Markov nature of the coupling and of the matching of the phonon wavelength to the 
interdot distance. 

PACS numbers: 71.38.-c,78.67.Hc,03.67.-a 
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Semiconductor quantum dots (QDs) are often consid- 
ered as candidate devices for a solid-state implementation 
of quantum information processing [lj, |2j, |3, La] • This in- 
terest is justified by the alleged long-lasting coherence of 
the interband polarization, which might be exploited for 
storing quantum information. Robustness against deco- 
herence in QDs is usually expected as a result of quan- 
tum confinement, implying discrete energy levels below 
the semiconductor bandgap with a very restricted phase 
space available for various scattering mechanisms. This 
simplified view was however recently questioned, after 
the observation of a very efficient decoherence mechanism 
due to strong electron-phonon coupling 0, 0, 0- The 
strongly localized polarization in an excited QD can in- 
duce virtual inelastic phonon emission and reabsorption 
processes which act as a phase-destroying mechanism. 
Theoretically, this stems from the exact solution of a two- 
level system coupled to a phonon bath - the so called in- 
dependent boson model In practice, the decoherence 
rate of the interband polarization is extremely fast in the 
first few picoseconds following excitation and vanishes af- 
terwards. To this corresponds a distinctive spectral fea- 
ture with a sharp zero-phonon line (ZPL) and phonon- 
assisted broad sidebands. Both time-resolved and spec- 
tral signatures of this mechanism have been characterized 
by several groups 0, H H H EH El E3 

According to the independent boson model, the ZPL 
relative to a single QD electron-hole level is not broad- 
ened. The measured ZPL linewidth, and the consequent 
dephasing time in the nanosecond range Q, are seem- 
ingly related to the radiative recombination process E3 • 
A recent theoretical analysis E3 has however suggested 
that virtual phonon-assisted transitions to excited lev- 
els might provide a dephasing of the ZPL, with a rate 
comparable to the radiative one. 

Recently, the system of two vertically stacked QDs, 
called QD molecules, has been experimentally studied 
lialMlllI- A QD molecule might in principle provide 
the minimal system of two coupled q-bits required for 
the implementation of quantum gates 0, Q • The elec- 
tron and hole wave functions in QD molecules can tun- 
nel through the thin barrier separating the QDs, giving 



rise to two bright (and two dark) exciton states fl5t Eif , 
whose energies are split by the effect of tunneling through 
the barrier and by Coulomb interaction. In a recent 
experimental study u sing coherent ultrafast four-wave- 
mixing spectroscopy |l4| . the decoherence time of the 
exciton interband polarization in QD molecules has been 
measured. It turns out that the ZPL linewidths are gen- 
erally larger than in the case of a single QD by almost 
one order of magnitude. Furthermore, the linewidths de- 
crease dramatically as the interdot distance is increased. 
Quantitatively, these larger linewidths can only partially 
be explained in terms of the larger volume of exciton 
states in a QD molecule, while the strong dependence on 
the interdot distance is still unexplained. 

In this letter we present a model of exciton optical 
response in a QD molecule, accounting for the strong 
exciton-acoustic-phonon interaction. The QD exciton 
states arc described within the effective mass scheme and 
the electron phonon coupling is treated within the sec- 
ond order Born approximation in the non-Markov limit. 
This limit is known to reproduce qualitatively - and to 
a good extent quantitatively - the results of the inde- 
pendent boson model |lfjj |. Linear optical spectra are 
computed at varying interdot distance. The measured 
behaviour of the ZPL linewidth is well reproduced. The 
larger linewidths are demonstrated to originate from the 
phonon-assisted coupling between phonon sidebands of 
one level and ZPL of a distinct level, lying at the same 
spectral position. This mechanism is essentially analo- 
gous to the excited-level scattering E3) but much more 
efficient due to the spectral overlap. The strong depen- 
dence on the distance is instead attributed to a wave- 
matching effect between the phonon wavelength and the 
interdot distance, which enhances the phonon-assisted 
scattering from bright to dark states. 

The dot molecule consists of two vertically stacked 
cylindrical QDs at mutual distance d, as sketched in 
Fig. n We assume throughout this work that the two 
QDs have identical shape and size, with radius R = 10 
nm and QD height h = 2 ran. Electron and hole states 
are described within two-band effective-mass approxima- 
tion yjj. As material parameters for InAs QDs, we 
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FIG. 1: Electron-hole pair energies (dashed) and exciton ener- 
gies (solid line) as a function of the interdot distance d. Inset: 
the geometry of the QD molecule. 



assume finite band offsets V e = 380 meV and V h = 
200 meV between the QD and the barrier material. We 
take m e — 0.067 mo, m^x — 0.34 too, and to^i = 
0.11 too for the electron, vertical and in-plane hole effec- 
tive mass respectively. For the effective-mass calculations 
we adopt cylindrical coordinates {p, z). In this geometry, 
when considering only QD-confined levels, a very good 
approximation consists in a separation of in-plane and 
z variables. Given the large energy separation of higher 
levels, compared to the Coulomb and deformation poten- 
tial interactions, we can restrict our calculations to the 
two lowest electron and hole levels. Without Coulomb 
interaction, the valence and conduction wave functions 
can therefore be written as ipj^(r) = fy(p)hjy(z), with 
energies Ejj, where 7 = v, c and j = 1, 2 for the symmet- 
ric and antisymmetric state respectively. In Fig. ^ both 
the bare electron-hole pair energies Ej c — E\ v (j, 1 = 1,2) 
as computed within the effective mass scheme, and the 
resulting exciton energies once Coulomb interaction (see 
below) is included, are plotted as a function of d. The 
energy splitting at short d is mainly due to conduction 
electron tunneling between the two QDs and vanishes, 
for electron-hole pair energies, at larger d. The Coulomb 
interaction produces both a global redshift and an ad- 
ditional splitting of the exciton energies, this latter per- 
sisting at large values of d. With the present parameters, 
the exciton transition energies measured in Ref. are 
well reproduced. More detailed models of the electronic 
states show that the detailed energy level structure 
deviates from this simple description, especially concern- 
ing the valence states. The phonon-assisted effect we 



are modeling, however, turns out to depend mainly on 
exciton-phonon coupling and, due to the lighter effective 
mass, almost entirely on the conduction-electron part of 
this coupling. The hole wave function affects only min- 
imally our results, finally justifying our adoption of a 
simple effective mass scheme. 

Within this scheme, the total Hamiltonian is H = 
H + H c + H$ + fig* where H = £,„ l^C^A 
is the bare electronic contribution in terms of the Fermi 
creation and annihilation operators. We include the 
electron-hole Coulomb interaction in our model because 
the phonon-mediated effects we are describing depend 
on the difference in binding energies between distinct ex- 
citon levels. Neglecting intraband terms, the Coulomb 
Hamiltonian H c has the usual expression 

4 = jE vzc\ c c} v c mv c lc , (1) 

ijlm 

with the Coulomb matrix element given by 

V lm= f dr edrhipi C (r e )iP* jv (r h ) 6 ^ B ipi c {v e )^ mv {v h ) , 

J | r e r M 

(2) 

being £b = 12.5 the background dielectric constant for 
In As. 

The two QDs are coupled through electron-acoustic 
phonon interaction via the deformation potential inter- 
action. The diagonal and non-diagonal parts of the 
electron-phonon interaction are 

H { ; h d) = EE M ^ « + a-0 <5JA . ( 4 ) 

where the matrix element for deformation potential cou- 
pling with acoustic phonons of dispersion w q = qs is given 
by 



= D lS j J^-t J dr^ 7 (r)^ 7 (r)e^- , (5) 

being s the sound velocity, p m the mass density, V the 
normalization volume and D 1 the deformation potential 

constants. The diagonal term describes the coupling 
of each isolated electron or hole state with the phonon 
modes 0, El El ■ The off-diagonal term accounts 
for the phonon-assisted scattering between different elec- 
tron or hole states. For our calculations, we take s = 4557 
m/s, p m = 5.67 g/cm 3 , D c = —13.6 eV and D v = —7.1 
eV Q. 

The independent boson model can be diagonalized ex- 
actly in the case of a single electron- hole pair state ||. 
Here however, because of the off-diagonal coupling H^f^ 
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FIG. 2: Comparison between the non-Markov (solid) and 
Markov (dot-dashed) computed exciton absorption spectrum 
at fixed temperature T = 10 K and interdot distance d — 
6 nm. 



and of the Coulomb interaction, the full diagonaliza- 
tion of the problem is a cumbersome task. Our ap- 
proach therefore consists in first diagonalizing the elec- 
tron Hamiltonian Hq + H c , giving rise to exciton states 
$j(r e ,r h ) = J2i m a l ™^i c (r e )ip mv (r h ), with energy Ej 
(j = 1, . . . , 4). The exciton-phonon Hamiltonians © and 
Q can be consistently rewritten in this basis in terms 
of exciton operators Bj,Bj. We then develop a density 
matrix formalism for the linear interband polarization, in 
the framework of the second-order Born approximation 
for the exciton-phonon interaction, without performing 
the Markov limit. It is known that this approach can 
quantitatively fairly well reproduce the exact result in 
the single QD case 0. For comparison, we will evalu- 
ate the Markov limit as well. The end equation for the 
four excitonic polarizations IF, (a;) = (Bj(u)Bj(t = 0)) 
reads 

[Ej -h(uj + ijj)] ILj(uj) - ^ Sjm(w)n m (o)) = gj , (6) 



where the exciton-phonon self energy is given by 
SjmH = Ek,* M Z kD *\ u ) M % being Ml k linear 



combinations of 



in the Coulomb basis, and 



Dq (lu) are the phonon propagators |8| (we assume a 
thermal phonon distribution at the lattice temperature). 
Here, h'jj is the radiative linewidth of the j-th excitonic 
state, computed as in Ref. 0, and 9j = \ J ^j(r, r)dr\ 
its oscillator strength. We numerically solve Eq. © for 
the exciton polarizations IFj(w). The optical spectrum 
is then related to the imaginary part of these quantities. 
A typical exciton spectrum for d — 6 nm is plotted in 
Fig. |21 Two of the four exciton wave functions are even 
with respect to z-inversion and the corresponding transi- 
tions are therefore optically active, giving rise to the two 
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FIG. 3: (a) Computed ZPL linewidths as a function of the 
interdot distance d at T = 10 K. Solid line: non-Markov. 
Dashed lines: Markov limit, (b) Same as in (a) but for T — 80 
K. 



peaks in the calculated spectrum fl5t H(| . Each peak is 
characterized by a narrow zero-phonon line (ZPL) and 
by phonon sidebands originating from the diagonal el- 
ements Yih(oj) of the exciton-phonon self-energy in lpj|l. 
as in the single QD case || 0|. The two other exci- 
ton states are dark. All four exciton states, however, are 
coupled to each other via the off-diagonal terms of the 
exciton-phonon self-energy in JHJ . The oscillations of the 
sidebands are an artefact of the Born approximation and 
are known to be absent in the exact solution of the inde- 
pendent boson model [lfl |. For comparison, we also plot 
the lorentzian spectrum computed in the Markov limit. 

Fig. Gla) shows the computed ZPL linewidths of 
the optically active excitons as a function of d, to- 
gether with the corresponding result in the Markov limit. 
The exciton-phonon coupling has two unexpected conse- 
quences. The first is a large broadening of both ZPLs 
at short QD distance d, which quickly decreases as d is 
increased. The same trend was measured by Borri et al, 
and left unexplained. From our calculation we can in- 
fer that this linewidth enhancement is due to transitions 
between bright and dark exciton states via phonon emis- 
sion or absorption. Let us consider only the conduction 
electron part of the exciton state, as the phonon cou- 
pling is much less effective for the valence band due to 
the heavy mass. The two conduction states have oppo- 
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site parity with respect to spatial z-inversion. In order 
for the integral in Eq. (J5J to be large, therefore, the 
phonon wavelength has to be an integer multiple of the 
interdot distance, thus restoring the even parity in the 
argument of the integral. If this condition is satisfied to- 
gether with energy conservation, the scattering process is 
enhanced. For d = 4 nm, the bright-dark energy splitting 
is 1.6 meV which, according to the phonon dispersion, 
corresponds to a phonon wavelength A ~ 12 nm, exactly 
three times the interdot distance. This wave-matching 
condition is quickly lost for larger d, resulting in the dra- 
matic linewidth decrease. This effect would be equally 
predicted by a standard Boltzmann dynamics and there- 
fore appears also in the Markov limit. The second conse- 
quence of exciton-phonon coupling observed in the data 
of Fig. inta) are the large values of the ZPL linewidths 
persisting in the non-Markov case even at large interdot 
distance. A similar behaviour was observed in the exper- 
iment by Borri et al . The computed linewidths must 
be compared to the Markov result which, for large d, gives 
much smaller values. The Markov values at large d are 
essentially given by the computed radiative linewidths, 
which are very different for the two brig ht exciton states 
due to their different symmetry HI EJ. The large mea- 
sured linewidths are therefore explained as being a con- 
sequence of the non-Markov nature of the coupling and 
of the coupled dynamics described by Eq. |JBJ. Due to 
the close energy spacing, each ZPL is spectrally super- 
posed to the tail of the phonon sideband relative to the 
other exciton levels, including the dark ones. Via the 
off-diagonal elements of the self-energy in ®, the ZPL 
is thus coupled to a continuum of levels, resulting in a 
spectral line broadening. This mechanism is analogous 
to the broadening of the ZPL due to virtual phonon- 
assisted transitions to excited QD levels, recently de- 
scribed by Muljarov et al E3, which can account for the 
ZPL broadening measured for sing le QDs @. Due to the 
close energy spacing of the four exciton levels, however, 
we predict that this is the most effective line-broadening 
mechanism for ZPLs in a QD molecule. At low temper- 
ature, phonon sidebands are strongly asymmetric, with 
a much larger high-energy side due to phonon emission 
as compared to phonon absorption probability. There- 
fore, the linewidth of the high-energy ZPL, according to 
the present mechanism, is expected to be larger than the 
low-energy one. This is seen in Fig. G^a), where the 
two ZPL linewidths differ by approximately an order of 
magnitude in the non-Markov case. The same difference 
was experimentally observed |l4| and attributed to the 
larger probability of relaxation from the upper to the 
lower level through phonon emission. At higher tem- 
perature, the phonon sidebands become more symmetric 
and, as a consequence, our mechanism predicts similar 
ZPL linewidths for the two spectral peaks, as shown in 
Fig. E2b) where the computed ZPL linewidths for T=80 
K are plotted. This further supports our interpretation. 



We expect however that, for this range of temperatures, 
coupling to other states via both acoustic and optical 
phonons come into play, thus requiring an extension of 
the present model. 

In conclusion, we modeled the optical response of exci- 
tons in a QD molecule, including the coupling to acoustic 
phonons beyond the Markov limit. Two new effects, ab- 
sent in the single QD case, dominate here the ZPL broad- 
ening. The first is the phonon scattering to dark states 
which is enhanced when the phonon wavelength matches 
the interdot distance. The second is the off-diagonal 
phonon-assisted coupling which, in the non-Markov limit, 
produces a spectral broadening of the ZPLs significantly 
larger than the radiative linewidth. All experimental 
findings can be explained in the the context of the present 
model. These broadening mechanisms provide, for QD 
molecules, a principial upper bound to the characteristic 
dephasing times, and represent therefore an even more 
restricting limitation, than in the single QD case, to the 
applications in quantum information technology. 
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